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> ' Abstract 

o 

(f-) ■ The driving terms in three-body theories of elastic and inelastic scattering 

O' 

of a charged particle off a bound state of two other charged particles contain 

ON 

the fully off-shell two-body Coulomb T-matrix describing the intermediate- 
state Coulomb scattering of the projectile with each of the charged target 
particles. Up to now the latter is usually replaced by the Coulomb poten 



tial, either when using the multiple-scattering approach or when solving 



three-body integral equations. General properties of the exact and the 
approximate on-shell driving terms are discussed, and the accuracy of this 
approximation is investigated numerically, both for atomic and nuclear pro- 
cesses including bound-state excitation, for energies below and above the 
corresponding three-body dissociation threshold, over the whole range of 
scattering angles. 
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I. INTRODUCTION 



Elastic and inelastic scattering of a particle off a two-particle cluster is conveniently 
formulated in terms of the exact three-body integral equations in momentum space (Fad- 
deev 1960) which, after suitable manipulations, can be given the practical form of mul- 
tichannel equations of the two-body Lippmann-Schwinger type (Alt et al 1967). If the 
impinging, and one or both of the target particles, are charged, modifications of this the- 
ory are required (Veselova 1971,1978, Alt 1978, Alt et al 1978,1980). Nevertheless, the 
basic, effective-two-body, multi-channel structure of the equations can be preserved. 

An alternative approach for calculating the (in-)elastic scattering amplitude, which is 
well-suited especially for high energies, is based on the multiple-scattering expansion of 
the corresponding three-body transition operator. Obviously, there exists a close corre- 
spondence between the terms occurring here and the various contributions to the effective 
potential of the integral equations approach. 

In applications of either formalism to charged-particle scattering, terms occur which 
contain the fully off-shell two-particle Coulomb T-matrix, describing the intermediate- 
state Coulomb scattering of various charged two-body subsystems. It is evident that the 
singularity structure of the latter in momentum space renders the calculation of such 
terms a rather difficult task. Hence, in numerical work (see, e. g., Alt et al 1985, Alt et al 
1994, and references therein) the Coulomb T-matrix is frequently replaced by its Born ap- 
proximation, namely the Coulomb potential, which drastically reduces the analytical and 
numerical effort required. But the reliability of such an approximation, in the following 
called 'Coulomb-Born approximation', is difficult to estimate. 

The first numerical studies of the quality of such a Coulomb-Born approximation 
appear to have been performed by Sinfailam and Chen (1972), who calculated the lowest- 
order multiple-scattering contributions to the elastic scattering amplitude for various 
atomic three-body processes, with the full two-body Coulomb T-matrix and with the 
latter replaced by the Coulomb potential. It turned out that, for most of the three-body 
systems studied, and for the energy and angular ranges considered, the Coulomb-Born 
approximation failed rather dramatically, the more in fact the lower the energy was. Of 
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course, if the latter is sufficiently high, both the exact and the approximate expression 
are expected to eventually yield the same answer; this was, indeed, found although at 
an energy larger than expected. This investigation, however, concerned only a limited 
variation of the scattering angle and the (altogether rather high) three-body energy. In 
particular, energies in the neighbourhood of the dissociation threshold, which represent 
critical tests of the reliability of the calculational procedures since the essential singu- 
larity of the two-body Coulomb T-matrix at zero subsystem energy plays a particularly 
important role there, were not considered at all. 

Theoretical and numerical studies of the lowest-order (off-shell) driving terms, which 
occur in the three-body approach, have been performed by van Haeringen (1979), and 
Kok et al (1979,1980,1981,1982). However, the investigations were restricted to negative 
total three-body energies, i.e., to energies below the composite-particle breakup threshold 
(and to the case that the masses of all three particles are equal, and the charges of the 
two particles involved in the intermediate-state Coulomb scattering are of equal sign). In 
some situations, which relate more to nuclear reactions with light nuclei, the Coulomb- 
Born approximation was found to be reasonably accurate. But in other cases, which are 
more typical for the situation prevalent in atomic physics, its lack of accuracy considerably 
diminishes its usefulness, similarly to the situation found at high energies. 

In this paper we address ourselves to a systematic investigation of the quality of the 
Coulomb-Born approximation, by comparing the exact with the corresponding approxi- 
mate on-shell driving terms for elastic and inelastic scattering in three-body systems. In 
particular, we study the dependence on several important parameters. Namely, energies 
are considered from the reaction threshold up to such high values that the Coulomb-Born 
approximation practically yields the exact result. Scattering angles are varied over the 
whole angle regime, with particular emphasis on very small angles which are favoured 
in higher-energy scattering. Furthermore, different masses of the two particles which ex- 
perience the intermediate-state Coulomb scattering, are allowed for: two light, one light 
and one heavy, and two heavy particles. Also the case that the magnitude of the charge 
of one of these two particles is larger than one is considered. And, finally, we also study 
the collisional excitation of the bound system. However, in the present investigation we 

3 



restrict ourselves to repulsive intermediate-state Coulomb interactions. 

In Sect. |l]we briefly sketch the three-body approach to elastic and inelastic scattering 
and, in particular, introduce the relevant driving terms. Some known results for the two- 
particle off-shell Coulomb amplitude which are used in the subsequent investigations are 
collected in Section [III]. In Sect. [TV], general properties of these elastic and inelastic 
driving terms are discussed, both in their exact form as well as if the intermediate-state 
Coulomb T-matrix is replaced by the Coulomb potential. Numerical tests of the accuracy 
of this approximation are presented for several typical atomic and nuclear three-body 
systems in Sect. [V], both for energies below and above the corresponding bound state 
dissociation thresholds. The results obtained are summarized in Sect. |VI| 

Natural units H — c — 1 are chosen. Furthermore, the conventional notation for 
two-body quantities: A a = A@ 7 , with a ^ [3 ^ 7 7^ a, is adopted. 
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II. THREE-PARTICLE MODEL OF ELASTIC AND INELASTIC SCATTERING 



Consider three particles with masses m v and charges e u , v = 1,2, 3. We are interested 
in the reaction a + (/?7) m — > a + {f3 r y) n) where particle a, having center-of-mass momentum 
q Q , impinges on the bound state of particles j3 and 7; the bound state wave function, 
belonging to the binding energy E am (with quantum numbers m), is denoted by \ip a m)- 
In the final state, the bound state quantum numbers are characterized by the index n, 
and the center-of-mass momentum of particle a is q' a . 

The corresponding (in-)elastic scattering amplitude T anjam (q' a , q a ; E) can be written 
as matrix element of an effective-two-body operator T anam (E + iO) between the plane 
wave states describing the free asymptotic motion of the clusters, 

(q' a ,q a ;E) = (q' a \T arijarn (E + iO)\q a ), (1) 

with the total three-body energy E being connected with the incoming and outgoing 
momenta by the on-shell condition 

E = ^- a +^ = ^- + E an . (2) 

Here, M a = m a (mp + m 7 )/(m Q + mp + m 7 ) is the a-channel reduced mass. 

According to the effective-two-body formulation of the three-body scattering theory 
(Alt et al 1967), the operators T antCem (z), together with the corresponding operators 
Tpn,am{z) appropriate for rearrangement scattering, are given as solutions of the following 
set of coupled Lippmann-Schwinger-type equations 

3 

T(3n,am{ z ) = Vpn,am{z) + Vfln/yr { z )Gor/,rs (^j^s.am (z) • (3) 

7=1 rs 

The sums over r and s extends at least over all bound states of subsystem 7, consisting of 
particles a and (3. Thus, the solution of (|3|) yields simultaneously the transition amplitudes 
for all two-fragment reactions. We note in parentheses that in the presence of Coulomb 
forces, the standard methods of integral equations theory can not be applied directly to 
(0) to obtain the physical reaction amplitudes, due to lack of compactness of their kernel 
(for a proper procedure in this case see Alt (1978), Alt et al (1978,1980)). 



In leading order in an iterative solution of eq. (|3|), the on-shell amplitudes describing 
elastic and inelastic processes are given by the on-shell 'diagonal' effective potential 

Tan,amiS\ a i Q.f*: 

= (c4KVan|f/4„(^ + iO)|Vam>|q«>. (4) 

The auxiliary three-body operators Up a satisfy the equations 

u' Pa {z) = 5 Pa G \z) + £ ^ 7 r;(^)Go(^)c/; Q (^). (5) 

7=1 

Here, Gol- 2 ) = ( z — H Q )~ l denotes the free three-body resolvent, and 5p a = 1 — 5p a the 
anti-Kronecker symbol. T'(z) is the sum of the two-particle off-shell Coulomb transition 
operator T^{z) plus a 'remainder short-range transition operator' in subsystem 7. For the 
definition of the effective free Green function Go-,-y,rs{z) and further details we refer to Alt 
(1978), Alt et al (1978). We only mention that the nondiagonal effective-potential parts 
Vf3n,am{ z ) with (3 7^ a, which are the driving terms relevant for rearrangement scattering, 
are defined similarly to @ in terms of the operators U'* a (z). They are, however, not 
considered in the present investigation. 

A Neumann series expansion of ([|) leads to the following representation of the 'diag- 
onal' on-shell effective potentials 

3 

Van, am 

(q' a , q a ;£) = 5 ia {q[ a \{ip an \T^ (E + iO)|^ am )|q a ) H . (6) 

7=1 

Explicitly shown are here only the leading Coulomb contributions due to (Coulomb) single- 
scattering in the intermediate state, which will be dealt with in the following investigation. 
The dots indicate single-scattering contributions from the remainder short-range transi- 
tion operator as well as higher-order rescattering terms. We note in parentheses that, if 
only two of the three particles are charged, and if the shorter-ranged interactions are rep- 
resented in purely separable form, the infinite series for V an , am (q' a , q a ; E) collapses exactly 
to just the terms written down explicitly in @. Furthermore, the terms @ would occur 
as the Coulomb single-rescattering contributions to the (in-)elastic scattering amplitude 
in a multiple-scattering approach. 

Let us introduce for the terms on the right-hand side of (|(J) the notation 
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KnMa,^,E) : = (q' a \^ an \T^(E + tO)\^ am )\q a ) 

= J (0Cn(p'J T^ P ' V p 7 ; E + iO - ^ am (Pa)- (7) 

The various subsystem momenta are defined as 

p Q = e ai I k + — q a I , = e a7 f k + — I , 

p 7 = e 7Q I q a + — k I , P 7 = e 7 « I q« + — k I . (8) 

Here, \i a = mpm 1 / (mp + m 7 ) is the reduced mass of the pair (Pj). For convenience, 
the antisymmetric symbol ep a = —e al 3, with e a p = +1 if (a, P) is a cyclic ordering of 
(1,2,3), is used. Moreover, the Coulomb T-matrix when read in the two-particle space is 
characterized by a hat, T c . The graphical representation of VJ^ m is given in Fig. 1. 

Similarly we define the quantities V^ m (q' Q ,q Q ) by the expression (0) but with the 
Coulomb T-matrix replaced by the Coulomb potential V? . They will be referred to 
as the 'Coulomb-Born approximation'. We mention that for simple types of bound state 
wave functions V^ m (q^,q Q ) can be calculated analytically (Lewis 1956). 
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III. THE TWO-PARTICLE COULOMB T-MATRIX 



In the literature there exist several equivalent expressions for the fully off-shell two- 
body Coulomb T-matrix in momentum space. A fairly complete account of results con- 
cerning two-body Coulomb scattering can be found in van Haeringen (1985). We will take 
the relevant formulae from van Haeringen (1979), and Kok and van Haeringen (1980) (see 
also van Haeringen and van Wageningen (1975), Chen and Chen (1972)), which were de- 
veloped there with regard to their applicability for the evaluation of expressions like the 
effective potential ([?]). Recall that we confine ourselves in the present investigation to case 
that the two particles participating in the intermediate-state rescattering have charges of 
equal sign. 

The following T-matrix representation is used, 

f 7 c (p', p; E, + zO) = V 7 C (A) jl + M| , ( 9 ) 

with 



^(A) = ^p, (10) 
I(y)=F irJi (y)-F tVi (y- 1 ). (11) 

Here, A = p' — p is the momentum transfer; E~ is the center-of-mass energy of the 



interacting pair of particles a and j3, and 7/ 7 = e a ep\J ^jlE^ the appropriate Coulomb 
parameter. Moreover, we have introduced the abbreviations 

x 2 1 (p /2 ~ 2/z 7 £ 7 )(p 2 -2 / i y £ 7 ) 
2/x 7 £ 7 A 2 

y = — 7> 13 

x — 1 

F ir]J (y) = 2 -Pi(i, "7 7 ; l + ivy, y), ( 14 ) 

with 2 Fi(a,b; c;y) being the hypergeometric function. Note that in writing down the 
formulae of this section we always assume E 7 ^ (formulae suitable for E^ in the vicinity 
of zero can be found in van Haeringen (1979)). For 3?z^ 7 > the integral representation 
for the hypergeometric function is particularly useful, 
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F irh (y)=i7 h £dtt*'- 1 {l-yt)- 1 , Mir] 7 >0. (15) 

It leads to the following expression for the Coulomb T-matrix 

{ r 1 t iv ~> 1 

TC(p>, p; E 1 + iO) = V°(A) h - Ai Vl J o dt u _ ^ 2 _ _ I , Mi Vl >0. (16) 

Furthermore, for < 0, which implies |y| < 1, one has the series representation for the 
function I(y): 



sinh 7rr] 7 7 ^— J n 2 + 77. 



We mention that off the energy shell, in the limit of vanishing momentum transfer, 
the following behaviour is easily derived: 

I(y) a^o 



x 



A{ci + c 2 lnA + c 3 A + o(A)}, (18) 



that is, the most singular part of the repulsive Coulomb T-matrix is given by the potential, 
for all (off-shell) subsystem energies E y . This leads for the ratio 1Z = T^(p', p; E 1 + 
iO)/V^ (A) to the well-known result 

\n\ A =° 1 + O(A) V E 1 . (19) 

Of course, on the energy shell one has \R\ — 1 for all scattering angles and energies. 

Kok and van Haeringen (1980) (see also van Haeringen and Kok (1984), van Haeringen 
(1985)) have derived useful general bounds on 71. Introduce the notation cos-# = p'-p/pp'. 
For instance, they have proved the following inequalities valid for a repulsive Coulomb 
interaction: 

0<TZ<1, V p,p', costi for £ 7 < 0, e a e p > 0, (20) 
0<|7£|<1, V p,p, costf for K, > 0, e a e >O, (21) 

which contain ( |I~9"D as special case. Illustrative graphical representations for the ratio 1Z 
are given in Kok and van Haeringen (1980). 
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IV. GENERAL PROPERTIES OF V^ nm AND V^ mn 

In this section we will discuss some general properties of the lowest-order contribu- 
tions to the elastic effective potential, graphically represented by Fig. 1, and investigate 
the effects which arise from the replacement of the Coulomb T-matrix by the Coulomb 
potential. 

Introducing the Coulomb T-matrix representation (|9|) into expression ([?[) one sees that 
the Coulomb potential can be taken out of the integral. Hence we can write 

V 7 T L(q^, q a , E) = Vf (A a )F^ m (q' a , q a] E), (22) 



where 



F^ m (qi,q a ;E) = J ^^n(p'J U + ^JVw(Pa), (23) 



with 



x 1 = 1 + 



pi - 2/i 7 (£ + i0- k 2 /2M 1 )} fp' 2 - 2// 7 (£ + iO - k 2 /2M, 



A a = e 7Q (q^ - q a ). (25) 
y is defined as before, cf. eq. (|T3|), and V!? by ([K]). Similarly, V^ m can be written as 

V 7 V l(qL, q«) = V^A^Ma, q«), (26) 



with 



/d^ k 
^C( k )im(k + fcA a /m„) (27) 

for m = n being called bound state form factor, and for m ^ n transition form factor. In 
the following investigation we will mainly be concerned with the ratio 

Ih^mVQa, <la, '■- yyc 7—, \ - F yc ( , \ ■ \ Z °) 

Let us state some general results for -R 7i „ m . 
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(i) For elastic scattering, both V^° mm and V^ m diverge for q„ — q Q — > because of (g) 
and ([H]); however, as eq. @ with flIB|) implies, their ratio tends in magnitude towards 
the value one, 

\Rr,mMa,^a,E)\ -> 1 for q^ q Q V £. (29) 



(ii) The inequalities (p0|) and (pTJ ) result in the following inequalities for the elastic 
ratio .R 7i oo; where the index denotes the ground state: 

0<i£ T)00 (c4q«;£)<l, for E<0, (30) 
< \R 7fi0 (q' a ,ci a ;E)\ < 1, for E > 0. (31) 

That is, for elastic scattering off a target in the ground state the Coulomb-Born approxi- 
mation always overestimates the exact driving term. Or, in other words, the error made 
by approximating in ([7p the two-body Coulomb T-matrix by the Coulomb potential is of 
known sign. Note that no analogous bounds result if either one or both bound state wave 
functions have nodes. 

(hi) Since for large two-body subsystem energies approaches the Born approxima- 
tion we expect for elastic and inelastic scattering 

R 7 , nm {<£ x ,q a ;E) E -^l. (32) 

However, it is obvious that for ( |32"D to hold the energy E must be higher than that 
for which on the two-body level we have T^(E y ) . For, in V^ m (q^, q a ; E) the 

Coulomb T-matrix enters for all 7-subsystem energies from E down to minus infinity, 
—00 < E 7 = E — k 2 /2M 1 < E. Thus a behaviour like (|3~2"1) can result only as a combined 
effect of T^'(E — /c 2 /2M 7 ) being approximately equal to VI over the whole range of 
momenta k for which the momentum-space bound state wave functions differ appreciably 
from zero. 

(iv) For excitation we deduce from ([27j) 



, lim FwJria,^) = for n^m (33) 
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because of the orthogonality of the bound state wave functions. Of course, on the energy 
shell (0) the momentum transfer can never vanish. However, in forward direction its mag- 
nitude, though remaining non-zero, can become small, thereby making R^ jnm (q' a , q a ; E) 
rather large. Furthermore, since the inequalities (|20| ) and fl2~I|) do not result in inequalities 
like ( |5DD or ( BT| ) for the ratio -R 7i „ m , values smaller and larger than one are possible. 



12 



V. NUMERICAL RESULTS 



When performing the integration over the magnitude of the momentum k of the non- 
interacting particle 7 in the exact driving term (|7|) we must, for a given three-body energy 
E > 0, treat the two regions 2M 7 £ < k 2 < 00 (region I) and k 2 < 2M^E (region II) which 
correspond to negative and positive subsystem energy f? 7 = E — fc 2 /2M 7 , differently. In 
region II, choosing 9ft 2 < 0, the representation ([Lj]) for the two-body Coulomb T-matrix 
can be used. In order to select an expression suitable for region I we first note that the 
Coulomb parameter r/ 7 occurring in depends on the integration variable k via 

r/ 7 = e a e /3 /i 7 /y / 2/i 7 (E + iO - k 2 /2M 1 ). (34) 

Thus, 7] 1 is imaginary. Since, as mentioned at the beginning, we restrict ourselves to that 
part of the elastic effective potential which contains Coulomb rescattering between the 
projectile and the target particle whose charge is of equal sign, i.e., e a ep > 0, in region I 
the condition 9ft irj^ > is always fulfilled so that the representation ([16]) is applicable. Of 
course, if expression ([7]) is calculated for a negative three-body energy E < 0, only region 
I is relevant. 

In order to show that the numerical quadratures can be performed reliably even for 
positive three-body energies (where the integration contour touches the essential singular- 
ity of the Coulomb T-matrix at zero subsystem energy), we present in Fig. 2 as a typical 
example the real and the imaginary part of (0), with all integrations except the one over 
the magnitude of k having been performed, as function of k. We choose the elastic reac- 
tion H(e, e)H to be discussed below, for a center-of-mass projectile energy of 100 eV and 
a scattering angle of 20 degree. Inspection reveals that the real part of this integrand and 
its first derivative are very smooth in the neighbourhood of that momentum for which the 
subsystem energy E — fc 2 /2M 7 = 0, in spite of the fact that, as discussed above, on either 
side of this point a completely different expression for the two-body Coulomb T-matrix is 
used. A similarly smooth behaviour is observed for the imaginary part. For other energies 
and scattering angles, and for other projectiles, qualitatively similar results are obtained 
(or course, if the projectile energy is smaller than the magnitude of the binding energy 
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of the incoming bound state, which corresponds to negative total three-body energy, the 
imaginary part of the integrand is identically zero). For comparison we also include the 



We furthermore mention that it was carefully checked that the driving term (|7|) behaved 
smoothly and yielded the same value when the ionisation threshold was approached from 
below and from above. 



For definiteness we consider only reactions where the bound states in the initial and in 
the final state are described by hydrogen-like wave functions. Let us begin by assuming 
the bound states to be in their ground states before and after the collision, i. e., we 
take m = n = {100} (with the usual notation {n£m} for the set of hydrogenic quantum 
numbers). Thus 



and the same expression for the outgoing bound state. Here, n 2 aQ = —2fi a E a0 = (/i a e/3e 7 ) 2 , 
E a o being the Coulomb ground state energy of the bound pair (fij). Introducing the wave 
functions ( |35"D into expressions ([?D and (|26f ) for the exact and the approximate driving 
terms, we can integrate analytically over the azimuthal angular integration variable. The 
remaining three-dimensional integrals have to be done numerically. It is to noted that for 
a few selected cases we have verified that our results coincide with those of Sinfailam and 
Chen (1972) within the accuracy with which the numbers can be extracted from their 
figures. 

In Fig. 3 we present the ratio \R eet00 (q' a , q a ; E)\ for elastic scattering of electrons off 
hydrogen atoms in their ground state, as function of the center-of-mass scattering angle 
■d (with cos $ = q a ■ q' a /q a Qa) an d of the center-of-mass projectile energy, starting at the 
elastic threshold. Here and in the following, the notation is such that the two particles, 



analogous (real) integrand of the Coulomb-Born approximation ( |2B| ) which, as is to be 
expected from is so much larger in magnitude than the real part of the exact expres- 
sion that the absolute value of the ratio is smaller than one, in accordance with fl3~ip. 



A. Atomic reactions 




(35) 
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which undergo Coulomb scattering in the intermediate state and which are denoted by 
the subsystem index 7 in fl28f ), are explicitly indicated. Furthermore, the ground state 
quantum numbers {100} in the initial and the final state are abbreviated by the index 
'0'. As is apparent, for each energy \R eet00 (q' a , q a ; E)\ starts in the forward direction at 
the value one, thus satisfying the condition (pup. For energies beyond 100 keV the ratio 
becomes practically equal to one again, in accordance with (|3^). Away from these regions 
this figure allows us to estimate the error made when using the Born approximation 
for the Coulomb T- matrix in the driving term (|7|). Inspection reveals that over a wide 
range of angles and energies the Coulomb-Born approximation dramatically overestimates 
the exact rescattering contribution. It is interesting to note that the minimum value of 
|i? ee) oo(qQ, q Q ; E)\ of approximately 0.27 is reached not in the backward direction but at 
an angle of about 60 degree (for a projectile energy of around 70 eV). 

The effect of increasing the mass of one of the particles undergoing intermediate- 
state Coulomb scattering can be inferred from Fig. 4. There we show the ratio 
|i? pe + i00 (q^, q a ; -E 1 )! for positron-hydrogen scattering, with the positron scattering off the 
proton. It is apparent that the picture is similar to the previous case, but minimum has 
deepened (to about 0.054 around 60 eV) and occurs in backward direction. 

The case of two heavy particles experiencing Coulomb scattering in the intermedi- 
ate state is presented in Fig. 5, which contains the ratio \R PP fi (c^ a , q a ; E)\ for proton- 
hydrogen scattering. Evidently, the mass effect is very pronounced, the Coulomb-Born 
approximation failing completely nearly everywhere, even at extremely small scattering 
angles for higher energies. The inadequacy of the Coulomb-Born approximation becomes 
still more pronounced if one of the heavy particles has a charge of magnitude greater than 
one. This can be inferred from \R pC <i+flo(<i' a , q^; E) \ for the elastic reaction C 6+ + H(ls), 
with C 6+ considered as a structureless particle, which is displayed in Fig. 6. In both cases 
the minimum (of the order 10~ 6 or less) occurs at 180 degree. 

In order to study the quality of the Coulomb-Born approximation for excitation we 
calculated the exact and the approximate driving terms (0) and (|26|) by assuming that 
the outgoing bound state is in a 2s-state, characterized by a wave function 
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VUP) = V>20o(p) = 64(2< )^ ( ^ + ^° )3 . (36) 

To simplify the notation, the bound state quantum numbers {200} will be abbreviated 
by the index '1'. We have already pointed out in Sect. [IV] that for excitation the mag- 



nitude of the ratio fl28|) is no longer bounded from above by the value one. In fact, the 
Coulomb-Born approximation becomes very small in forward direction, thereby giving rise 
to large values of the ratio. This situation is exemplified in Fig. 7 where we have plotted 
|-Ree,oi(qa, Qa! E) | for excitation of hydrogen atoms from the Is- to the 2s-state by electron 
impact. It is seen to take on values larger and smaller than one. The value one is again 
reached at higher energies similar to those for elastic scattering shown in Fig. 3. Fur- 
thermore, although barely visible in the figure, \R eet oi(q' a , q a ; E)\ develops a minimum of 
0.27 at 45 eV projectile energy and 90 degree. Interestingly, this minimum value coincides 
with that for the elastic ratio and, as was the case there, occurs at non-backward angles. 
Since, when approaching the forward direction, |i? eej01 (q^,, q Q ; E)\ increases sharply, this 
figure has been cut off at a scattering angle of 1 degree. It is clear that the larger the 
energy is, the relatively smaller values of the momentum transfer are accessible and thus 
the smaller the Coulomb-Born approximation becomes (cf. eq. (|33|)). Of course, for suffi- 
ciently high energy the near- vanishing of V^ < g 1 (q^,, q a ; E) in the forward direction has to 
be counterbalanced by a corresponding near-vanishing of V^ C 01 (q^,, q a ] E), in order that 
their ratio eventually approaches the value one. We only mention that a similar picture 
is obtained for the ratio \R pp fii(q' a , q a ; E)\ for excitation of hydrogen atoms from their 
ground to the same excited state by proton impact. 

Since due to these kinematic effects the ratio (|28| ) is no longer an appropriate quantity 
to visualize the quality of the Coulomb-Born approximation (except for scattering angles 
larger than, say, 90 degree), we plot in Fig. 8 only the absolute value of the exact 
driving term Vj^ 01 (q' a , q a ; E) for excitation of hydrogen atoms by electron impact, as 
function of the scattering angle and the projectile energy. It is a very smooth quantity, 
however ranging over 14 orders of magnitude in the range of parameters considered. The 
mass effect is again illustrated by considering the exact driving term Vj p 01 (q' a , q a ; E) 
for excitation of hydrogen atoms by proton impact, displayed in Fig. 9. Not only is the 
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range of accessible values smaller than in the previous case; but there appears also to be 
more structure in the amplitude. We point out that, at small scattering angles, the exact 
driving terms for both reactions increase in magnitude with increasing the energy, and 
decrease only at rather high energies; the energy where the maximum is reached lies the 
higher the smaller the angle is. 



In order to investigate whether these rather large deviations from the value one of 
|-R 7 ,oo| do hold also for other systems, we considered two examples of elastic deuteron- 
nucleus scattering, i.e., one of the three particles is a neutron. The nuclear interactions 
are assumed to be given in separable form. Recall that under such conditions, the expres- 
sion (0), with 7 denoting the proton-nucleus subsystem, represents the exact 'diagonal' 
contribution to the on-shell effective potential. The deuteron bound state wave function 
used here is the same as that considered, e. g., by Alt et al (1985) in their calculation of 
elastic proton-deuteron scattering observables. That is, the deuteron is assumed to be in 
an pure S state, with a momentum space wave function 



Here, = y —m N E d , with E d = —2.226 MeV being the deuteron binding energy and 
tun the nucleon mass; N is a normalisation constant, and (3 a parameter fitted to the low 
energy neutron-proton scattering parameters in the 3 Si-state. 

First we look at proton-deuteron scattering. With this bound state wave function 
we have calculated the ratio \R PP fl (q' a , q a ; E)\ with proton-proton intermediate-state 
Coulomb scattering, again as function of the center-of-mass scattering angle and proton 
bombarding energy. Note that the threshold for the deuteron breakup is at a center-of- 
mass proton kinetic energy of 2.226 MeV. The results are shown in Fig. 10. Clearly, 
\R PP fio(q' a , q a ; E)\ satisfies the general constraints (P9"D, fl3"U|), (|3"I|) and (^5J). But, as is 
apparent, it differs only very little from the value one, the maximal deviation of 12 per- 
cent occurring at energies between 6 and 7 MeV, in backward direction. The conclusion 



B. Nuclear reactions 




(37) 
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is that for such a reaction the use of the approximation — > VP in fl7|) appears to be 
well justified. 

Similarly to the atomic case, the Coulomb-Born approximation quickly ceases to be 
acceptable if the charge of one of the charged particles is increased. This is illustrated in 
Fig. 11, where we show the ratio \R p c G + floi^'ai Qai E)\ for the reaction of C 6+ colliding 
with deuterons (again, C 6+ is considered a structureless nucleus). Here, the ratio drops 
to 0.47 (at 12 MeV and a scattering angle of 67 degree). 
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VI. DISCUSSION 



We have investigated the quality of the replacement of the two-particle Coulomb am- 
plitude by the Coulomb potential in the driving term for elastic and inelastic scattering, 
both for atomic and nuclear three-body reactions. We have found the following interesting 
results. 

(i) For the atomic reactions studied, the Coulomb-Born approximation quite generally 
turns out to be completely unsatisfactory. Its failure becomes the more striking the heav- 
ier the masses of the particles involved in the intermediate-state Coulomb scattering are, 
and - not surprisingly - the higher their charge is. 

(ii) In the nuclear cases considered the situation is much more favourable. At least for 
proton-deuteron scattering the Coulomb-Born approximation, which on the energy shell 
has been shown to be generally accurate to better than 10 per cent, is fully justified (this 
is particularly so, since the effective-potential part (|7|) itself represents only a correction 
to the dominant driving term, which is induced by the nuclear interaction; in fact, its 
Coulomb-Born approximation fl26|) is known to contribute something of the order of 10 
per cent to the final proton-deuteron observables, cf. Alt et al 1985). Of course, also here 
its quality is diminished if the charge of one of the particles is increased, but the effect is 
not so dramatic as for the analogous atomic system. 

(iii) For atomic as well as nuclear elastic scattering processes it has been found that the 
exact driving term (|7|) becomes equal to the approximate one, eq. fl26|), not only in forward 
direction and for higher energies, but also if the projectile energy goes to zero or, in other 
words, if the elastic threshold is approached, for arbitrary scattering angle. We have not 
yet succeeded in developing a convincing argument for this unexpected fact. A first sug- 
gestion, however, is based on the known invariance of the ratio T!?(p', p; E^+iO) /V^(p', p) 
under the transformations p/k — > k/p, p'/k — > k/p', with k = \J 2m 7 i? 7 (cf, e.g., Kok 
and van Haeringen 1980), which correlate high-energy (k — > oo) with low-energy proper- 
ties (k — > 0). This speculation is supported by the fact that the value one of the ratio 
( [28] ) is reached for the larger and the smaller projectile energies, respectively, the larger 
the scattering angle becomes. 
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(iv) The minimum of the elastic ratio (]28f) is in all cases attained for energies above the 
dissociation threshold, frequently but not always in the backward direction. 

(v) For excitation the Coulomb-Born approximation is of similar quality as for the corre- 
sponding elastic reaction, except for the additional failure in the near-forward direction. 
Moreover, in this case the ratio does not approach the value one at the reaction 
threshold. 

(vi) In all reactions considered the ratio (|28|), and hence also the exact driving term ([?D, 
behave smoothly when the bound state dissociation threshold is crossed. 

(vii) In general, the imaginary part of the exact driving term (|7]) is much smaller than 
its real part (the Coulomb-Born approximation ( ^6|) is real everywhere), except in the 
region of parameters where the ratio approaches its minimum value. There they become 
comparable in magnitude. 

From the numerical point of view it may be of interest to note that the use of the 
Coulomb T-matrix in the calculation at least of the on-shell driving term (|7]) for parti- 
cles with charges of equal sign, does not present any problem with modern computers. 
Whether this holds true also for the calculation of similar expressions but with oppositely 
charged particles, or for the corresponding off-shell quantities (which would be needed as 
input in the three-body integral equations approach) remains to be seen. 
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Figure Captions 

Fig. 1. Graphical representation of the elastic and inelastic driving term V^ C nm . 
Fig. 2. The integrand of Vj 00 (q' a , q a ; E) (cf. eq. (g)) for e + H(ls) -> e + #(ls) as 
a function of the magnitude of the momentum k of the proton (the integrations over all 
other variables are already performed) at incident cm. energy of 100 eV and at cm. 
scattering angle 20 deg. Upper curve is the real part, and lower curve the imaginary part 
of the integrand. Also shown is the analogous integrand for the Coulomb-Born approxi- 
mation V^Q (q^, q a ; E), reduced by a factor of 5. The arrow indicates the point where 
the ee-subsystem energy is zero. The integrand and k are in a.u. 

Fig. 3. The ratio \R eet00 (q' a , q a ; E)\ as a function of cm. incident energy and cm. scat- 
tering angle for elastic electron + H(ls) scattering. 

Fig. 4. The ratio \R pe + fi (q' a , q a ; E)\ as a function of cm. incident energy and cm. 
scattering angle for elastic positron + H(ls) scattering. 

Fig. 5. The ratio \R PP flo(q' a , q Q ; E)\ as a function of cm. incident energy and cm. scat- 
tering angle for elastic proton + H(ls) scattering. 

Fig. 6. The ratio |-R p c 6 +,oo(qa5 q^i E)\ as a function of cm. incident energy and cm. 
scattering angle for elastic scattering of C 6+ off H(ls). 

Fig. 7. The ratio \R ee fii(q' a , q a ; E)\ as a function of cm. incident energy and cm. scat- 
tering angle for the inelastic reaction e + H(ls) — > e + H(2s). 

Fig. 8. The amplitude Vjf 01 as a function of cm. incident energy and cm. scattering 
angle for the inelastic reaction e + H(ls) —>■ e + H(2s). 

Fig. 9. The amplitude V pp01 as a function of cm. incident energy and cm. scattering 
angle for the inelastic reaction p + H(ls) —>■ p + H(2s). 

Fig. 10. The ratio \R PP fio(q' a , q a ] E)\ as a function of cm. incident energy and cm. 
scattering angle for elastic proton-deuteron scattering. 

Fig. 11. The ratio |-R p c 6 +,oo(qa> q^j E)\ as a function of cm. incident energy and cm. 
scattering angle for elastic C 6+ — d scattering. 
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